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ABSTRACT 

We present new mid-infrared (5 — 35/im) and ultraviolet (1539 - 2316 A) observations of the inter- 
acting galaxy system Arp 143 (NGC 2444/2445) from the Spitzer Space Telescope and GALEX. In 
this system, the central nucleus of NGC 2445 is surrounded by knots of massive star-formation in a 
ring-like structure. We find unusually strong emission from warm H 2 associated with an expanding 
shock wave between the nucleus and the western knots. At this ridge, the flux ratio between H 2 and 
PAH emission is nearly ten times higher than in the nucleus. Arp 143 is one of the most extreme 
cases known in that regard. From our multi-wavelength data we derive a narrow age range of the 
star-forming knots between 2 Myr and 7.5 Myr, suggesting that the ring of knots was formed almost 
simultaneously in response to the shock wave traced by the H 2 emission. However, the knots can be 
further subdivided in two age groups: those with an age of 2-4 Myr (knots A, C, E, and F), which are 
associated with 8/im emission from PAHs, and those with an age of 7-8 Myr (knots D and G), which 
show little or no 8/im emission shells surrounding them. We attribute this finding to an ageing effect 
of the massive clusters which, after about 6 Myr, no longer excite the PAHs surrounding the knots. 
Subject headings: galaxies: interactions — galaxies: starbursts — galaxies:ISM — galaxies: individual 
(NGC 2445) 



1. INTRODUCTION 

The direct study of warm molecular hydrogen emis- 
sion from galaxies has been improving considerably 
through the advent of more sensitive space-bourne 
spec trometers on the Infrar e d Space Ob s ervato ry (e. 
g. iRigopoulou et all (|2002t ): iLutz et all ([2003])) and 




Fig. 1. — Overlay of the SL (sparse) and LL (complete) maps on 
a V-band image of Arp 143. North is up. 
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more recently on the Spitzer Space Telescope 5 . With 
the Infrared Spectrograph (IRS 6 ) (|Houck et all l200l 
on Spitzer (|Werner et alj [2004) it has been possible 
to improve upon systematic studies of the mid-IR 
pure-rotational lines of molecular hydrogen in a va- 
riety of environmen ts, from nearby normal galaxies 
(|Roussel et al.l [20071 ) . to galaxies in more extreme en- 
vironments, like g r oups and clusters (lOgle et alJ 120071 : 
lEgami et all [20061 iFerland et alJ [2003) and ULIRGS 
(|Armusl 1200a iHigdon et all l2005h (see also a review on 
molecules by lOmontl (|2007h ). 

Recently unusually strong mid-IR H 2 lines have been 
discovere d in the giant shock-wa ve structure in Stephan's 
Quintet (jAppleton et all 120061 ) , and in over a doze n 
low-luminosity radio galaxies (jOgle et alJ I2007L l2008h . 
where H 2 are often the strongest in mid-infrared wave- 
lengths. Strong H 2 equivalent widths were also found 
near t he "overlap" regi on in the Antenna interacting 
galaxy (jHaas et al.ll2005| ) . These systems have extremely 
large H 2 equivalent widths and line luminosities (L# 2 > 
10 41-42 ergs/s), and relatively low star formation rates. 
Another defining characteristic is the large ratio of H 2 
luminosity to PAH line strength, and the large (0.1 to 
30%) fraction of H 2 luminosity to the total bolometric lu- 
minosity of the objects. Even larger H 2 line luminosities 
have been found associated with several massiv e galax- 
ies in X-ray clusters (e. g. lEgami et al.l (|2006h ). This 
new class of powerful H 2 emitting galaxy appears to be 
shock excited, and models of Stephan's Quintet, where 

5 This work is based, in part, on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract 
with NASA. Support for this work was provided by NASA through 
a G02 award issued by JPL/Caltech. 

6 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 
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a large-scale shock is strongly implicated suggest that a 
significant fraction of the bulk kinetic energy in the shock 
must be funnelled into the H2 line in order to explain the 
results (|Boulanger et al1l2008t iGuillard et all 120081) . In 
this context, it is of considerable interest to find other 
examples of the same kind of emission in nearby sys- 
tems. This paper describes strong H2 emission from the 
strongly interacting galaxy pair Arp 143. 

Arp 143 is an interacting pair of galaxies (NGC 
2444/2445) with many of the necessary ingredients for 
a study of shocks generated by galaxy collisions. The 
southern gas-rich component of the pair, NGC 2445, is 
notable for the ring of stellar clusters and HII regions, 
whereas the SO companion is devoid of activity. It is pos- 
sible that NGC 2445 shares some similarities with colli- 
sional ring galaxies (like the Cartwheel ring), because the 
powerful star- forming knots lie along a cres cent-shaped 
wave of HI emission (jAppleton et al J 119921 ). and there 
is kinematic evide nce that the struct ure is expanding 
through the disk (iHigdon et alj 119971). Early unpub- 
lished spectroscopy by (jJeske et al.|1986l) and optical and 
near-IR photometry of the knots (jAppleton et al.lll992l ) 
suggest the star clusters are quite young <30-60Myrs, 
much younger than the dynamical age of the expand- 
ing HI wave-thus providing evidence that they a r e trig - 
gered by a collision with NGC 2444. Ueske et al] (119861 ) 
estimated sub-solar (LMC-like) metallicity in the disk 
knots in the range 12+log[0/H] = 8.56-8.77, similar to 
other known collisional ring galaxies (see Bransford et 
al. 1996). The nucleus has mildly super-solar values of 
12+log[0/H] = 9.18-similar to nuclear starbursts. 

However, this simple picture is complicated by the 
unsolved mystery of t he large 150kpc-long HI plume 
( Appl eton et al.lll987al ) which extends just to the north 
of NGC 2444. The existence of such a long plume 
implies that the two galaxies involved in the colli- 
sion have had a prev i ous major tidal encounter in the 
past dHibbard fc Yunlfl99l Iffibbard eTaT][2000T) . Per- 
haps the simplest interpretation of the system is that 
NGC 2444 may have experienced an initial close pas- 
sage with NGC 2445 ~100 Myrs ago, stripping HI 
into the plume from NGC 2444. After an initial non- 
pentrating encounter, the main stellar body of NGC 
2444 has returned to collide recently with the disk of 
NG C 2445 in a manner similar to co llisional ring galax- 
ies ([Appleton fc Struck-Mareelill 19961) . 

In this paper we will present new multi-wavelength 
imaging and spectroscopic data of Arp 143, another ex- 
ample of a galaxy with very strong Mid-IR H2 emission. 
Our goal is to connect the recent dynamical events in 
Arp 143 with the shock and star formation history of 
the galaxy. The paper will use Spitzer, GALEX 7 , and 
new ground based optical and near-IR images as well as 
mid-IR spectroscopy to explore the physical conditions in 
Arp 143. This will include the study of the dust and HII 
region properties, ages and luminosities of star formation 
regions and the excitation conditions of the warm H2 gas. 
After a description of the images and spectra in §2, we 
will present in §3 the multiwavelength photometry and 
the spectral analysis. In §4, the discussion of the results 

7 GALEX (Galaxy Evolution Explorer) is a NASA small ex- 
plorer launched in 2003 April. We gratefully acknowledge NASA's 
support under Guest Investigator program no. 45. 



will follow two main directions: the analysis of the shock 
region using mainly results from the spectral analysis, 
and the star formation history of the knots based on an 
interpretation of the UV-infrared SEDs of each knot. We 
will assume a distance to Arp 143 of 56.7 Mpc based on 
its corrected Virgo-centric velocity of 4142 km/s and an 
assumed Ho=73 km/s/Mpc. At this distance, 1 arcsecs 
corresponds to 275pc. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spectra 

Th e spectra were take n on October 24, 2004 using the 
IRS (jHouck et all 120041 ) on board of the Spitzer Space 
Telescope. All spectral data were processed using IRS 
pipeline version S15. Fig. Q] shows an overlay of the IRS 
low- and high-resolution slits on a V-band image. We 
have a completely sampled map in the Long-Low (LL; 
15 - 37 fim) spectral mapping mode, consisting of 17 
pointings (each with a 32 s exposure): a total LL expo- 
sure time of 544 s. A sparsely-sampled map was obtained 
with Short-Low (SL; 5-14 fxm) which did not fully sam- 
ple the galaxy as shown in the figure. It consisted of 
30 pointings with 61 s exposures each: a total SL ex- 
posure time 1891 s. Since the LL and SL observations 
are composed of two sub-slits, each observed separately, 
convenient off-source "background" exposures were au- 
tomatically obtained during the observations. These ob- 
servations were selected to be free of line emission and 
were subtracted from each on-source exposure. 

High-resolution spectra of the nucleus were taken in 
two nods with two 31 s exposures each for short and long 
wavelengths (SH and LH), resulting in a total integration 
time of 248 s. For each module we averaged the expo- 
sures for both nod positions. After the removal of bad 
pixels, we extracted the spect ra for each module, using 
SMART 8 (iHigdon et al.ll2004D . a tool for the extraction 
and analysis of Spitzer-IRS spectra, written in IDL. For 
the high-resolution spectra we did not subtract a back- 
ground since there was no suitable "sky" spectrum taken. 
These observations were made early-on in the Spitzer 
mission before it was realized that obtaining a "sky" ob- 
servation away for the target was useful for rogue-pixel 
removal (especially LH). Therefore, our SH and LH ob- 
servations contain low-level zodiacal light in addition to 
the continuum from the galaxy. This does not affect the 
measurements of line fluxes (see later) but potentially 
can affect the measurement of the absolute continuum. 
The source was sufficiently bright that rogue-pixel re- 
moval was not a large issue in the case. 

Extraction of all the in-targe t low-resolution spe ctra 
was performed using CUBISM (jSmith et al.ll2007bl ). an 
IDL tool for the construction and analysis of spatially 
resolved spectral cubes using Spitzer/IRS spectra. Bad 
pixels in the BCD images were manually flagged in each 
cube pixel, and then automatically discarded when re- 
building the cube. 

2.2. Images 

Fig. [5] shows images of Arp 143 at 12 different wave- 
lengths: far- (A 154nm) and near-UV (A 232nm) contin- 
uum; optical B- V-, and R-band continuum; H Q emission; 

8 SMART was developed by the IRS Team at Cornell University 
and is available through the Spitzer Science Center at Caltech. 
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Fig. 2.— Images of Arp 143 in 12 different bands: GALEX FUV, GALEX NUV, B, V, R, Ha, K, IRAC 3.6iim, IRAC 4.5/im, IRAC 
5.8/^m, IRAC 8/jm, and IRS LL map at 24 (im. The different apertures used for photometry are overlayed on the FUV image. 



K-band, the mid-infrared 3.6/im and 4.5^m bands, sen- 
sitive to the emission of the red stellar population; the 
mid-infrared 5.8/im and 8.0/im bands, specially sensitive 
to PAH and dust emission; and a 24^tm continuum image, 
built from a IRS LL spectral map. In the FUV image we 
show the apertures chosen for photometry, matching the 
main star forming knots visible on FUV, marked from 
A to G, and the nucleus. Each image was obtained as 
follows. 



2.2.1. FUV and NUV 

Ultraviolet images of Arp 143 were obtained on Jan- 
uary 28, 2005 (FUV) and December 26, 2006 (NUV) us- 
ing the Galaxy Ev olution Explorer (GALEX) satellite 
([Martin et al.ll2005[ ). The galaxy was imaged in the FUV 
and NUV bands, covering the wavelengths 1344 — 1786A 
and 1771 - 283lA. The total integration time was 3040 
seconds for the FUV image and 219 seconds for the NUV. 
GALEX uses two 65 mm diameter, microchannel plate 
detectors, producing circular images of the sky with 1.2° 
diameter at 5" resolution. The GALEX image was re- 
duced and calibrated through the GALEX pipeline. The 
fluxes were converted to magnitudes on the AB system 
([Okeill990h . The data was constructed with version 6 of 
the GALEX pipeline. We coadded two All Sky imaging 
visits to obtain the 219 seconds of exposure in the NUV 
image. 



2.2.2. Optical and Near-Infrared 

We obtained optical images on February 3, 2003 using 
the Palomar 60" telescope. The observations were made 
with a FOV of 12'. 9 x 12'. 9 (2048 x 2048 pixel) CCD 
imaging system with an RCA chip. The pixel scale is 
0".38/pixel The data were taken in three filters, John- 
son B, Johnson V, and Johnson R, and exposure times 
of 500 seconds for each filter. The night was moonless 
and transparent and the seeing was 0".7 — 1".0. We 
used the star RU 152 for calibration. The Ha image 
was published in I Romano et al.l ([2008( 1 and was taken 
in the Guillcrmo-Haro 2.1 m Telescope using a narrow- 
band filter centered at 66351 (FWHM ~ 97 A), which 
also contains the [NII]A6583 line. We corrected the Ha 
flux for [Nil] by assuming Ha/ [Nil] ~ 3, typical of HII 
regions (lOsterbrocld[l98ai . The near-infrared images 
were obtained with good (< 1") seeing through a K s 
(A = 2.15/im) filter using the WIRC camera on the Palo- 
mar 200" telescope in November 30 2004. WIRC is a 
2048 x 2048 pixel wide-field 8.5 x 8.5' HgCdTe camera 
with a pixel scale of 0.25"/pix operated by Caltech. 

2.2.3. Mid-Infrared 

Arp 143 was imaged with t he Infrared Array Cam- 
era (IRAC, iFazio et al.1 ([2004( 1 ) on Spitzer at 3.6, 4.5, 
5.8, and 8.0 ^m on March 28, 2005. The IRAC detec- 
tors consist of two 256 x 256 square pixel arrays with 
a pixel size of 1"22 resulting in a total field of view of 
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Fig. 3. — Low resolution spectra of six regions coinciding with 
the selected knots. 
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Fig. 4. — Low resolution spectrum of the gas ring in NGC 2445. 

5". 2 x 5". 2. For each of the four channels, 47 12 s ex- 
posures were taken, with a total integration time of 564 
seconds for each channel. The data were reduced us- 
ing standard procedures (i.e., dark-current subtraction, 
cosmic ray removal, non-linearity correction, flat-fielding 
and mosaicing) using pipeline version 14.0 of the Spitzer 
Science Center. 
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Fig. 5. — Full high resolution spectrum of the nucleus of NGC 
2445. The SH part was scaled up by 1.8 to match the LH spectrum. 



3. ANALYSIS 

In Fig. [3] are five low-resolution spectra of 21" x 21" 
regions coinciding with knots B, C, E, F, and the nu- 
cleus. In these spectra we can see broad features at- 
tributed to Polycyclic Aromatic Hydrocarbons (PAHs), 
such as 6.2^m, 7.7/xm, 8.6/mr, 11.3/xm, 12.6/im features. 
Also noticeable is the distinctive PAH 17/im complex. 
In the regions E and F, the H2 lines at 9.7/im, 12.2/im, 
17/xm, and 28.2^m are especially strong. The ionic lines 
observed are [Aril] at 6.9^m, [Nell] at 12.8/xm, [Nelll] 
at 15.6/im, [SIII] at 18.8/im and 33.5/im, and [Sill] at 
34.8/im. The spectrum of knot D is extremely noisy, and 
therefore it is not shown in Fig. [3l This means that there 
is no continuum flux detected at this wavelength. How- 
ever, the PAH features are detected above the sensitivity 
level. The peak of the 6.2/im PAH feature, for exam- 
ple, is at 4 mjy, well above the 3cr sensitivity of the SL 
slit at this wavelength, which is 0.6 mjy. In Fig. [4] is 
a low resolution spectrum from the region wher e a rin g 
of HI emission was observed by lAppleton et aTl (1992); 
iHigdon et alj (|1997l ). which is delimited in the B-band 
image in Fig. O This spectrum is very similar to the 
spectra of knots E and F, but it has some differences. 
It is even more dominated by distinctive H2 lines, espe- 
cially the H2 S(l) line at 17/itm, although PAH features 
can also be seen, especially at 11.3/im. The ionic lines 
are not so prominent, as the spectrum was taken in a 
region not dominated by star forming clusters. 

3.1. H2 lines 

Using the pure rotational H 2 lines we can probe the 
physical conditions of the warm molecular hydrogen in 
the star forming knots. The pure rotational lines origi- 
nate from the warm (T > 100 K) gas. We measured the 
H2 rotational line fluxes using SMART, and derived their 
excitation temperatures, column densities, and masses of 
their temperature components. Fig. [6] shows the excita- 
tion diagrams of the detected H2 transitions and upper 
limits for the knots E, F, G, the nucleus, and the ring. 
Each excitation diagram consist on a plot of the natural 
logarithm of the column density N divided by the statis- 
tical weight g in the upper level of each transition against 
the upper level excitation temperature T ex . The column 
density follows from the Boltzmann equation: 
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Fig. 6. — H2 excitation diagrams for knots E, F, G, the nucleus, 
(T), column densities N(H2) and molecular gas masses M(H2) in u 
are upper limits and are not included on the fit. 

N Z(T ex ) F \ T ex ) K ' 

where N(i) is the molecular column density of the ith 
transition, N is the total column density of H 2 , g(i) is 
the statistical weight for the ith transition, and Z(T ex ) is 
the partition function at the excitation temperature T ex . 
The values of g for odd and even transitions are differ- 
ent because of the ortho and para transitions. The H2 
line fluxes, and the molecular gas temperatures, column 
densities and masses are listed in Table Q] 

The excitation temperature of the line-emitting gas 
is the reciprocal of the slope of the excitation diagram, 
and corresponds to the kinetic temperature in local ther- 
modynamic equilibrium (LTE). However, the nonlinear 
decline of log (N/g) with upper-level energy, commonly 
seen in shocks within the Galaxy, as well as in exter- 



the gas ring . The plots are labelled with the fitting temperatures 
of K, cm -2 , and Mq respectively. The points that are not labeled 

nal galaxies (jLutz et alJ 1200 3: Rigopoulou et al.ll2002f ). 
is an indication that no single-temperature LTE model 
fits these data. 

The lowest temperatures are measured in knots E and 
G, using S(0) and S(l). The temperatures for knots F 
and H were derived from higher excitation transitions, 
S(l) and S(2), as S(0) was not detected. The measured 
temperatures are 137 K for knot E, 280 K for knot F, 
178 K for knot G, and 309 K for the nucleus. We can 
see that the molecular gas is hotter in knot F (280 K) 
and in the nucleus (309 K), but denser in knot E, where 
most of the warm molecular gas is concentrated. It is 
important to note that the temperatures for knots E and 
G were derived using only the S(0) and S(l) transitions, 
whereas higher transitions were used in other knots re- 
sulting in a higher average temperature of the gas. The 
derived excitation temperature is also dependent on the 
ortho-to-para ratio for H2. During the fit to a multi- 
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temperature model through the preferred data points, 
we constrain the warm gas by the S(0)/S(2) [para] ratio 
and allowing the higher order S(3) and S(5) [ortho] tran- 
sitions to provide a very rough guide to the temperature 
of a hotter component, resulting in ortho-to-para ratios 
of 2.2 for knot E, 2.7 for knot G, and - 3 for knot F, 
the nucl eus, and the ridge. These ratios agree with LTE 
models ([Burton et al.| [l992). The various masses are de- 
rived multiplying the column density by the physical area 
of the aperture corresponding to the knot. 

3.2. PAH features 

We used PAHFIT ([Smith et al.ll2007f) to measure the 
PAH emission features for each of the knots. Some PAH 
features are decomposed by PAHFIT into gaussian com- 
ponents, such as the 7.7/im and 11.3/im features, and 
merged for the calculation of their total flux and EWs. 
The results of the PAH feature fits are listed in Table (2[ 
Most of the PAH features lay on the SL part of the spec- 
trum, meaning that the sparse map only has a 2 pixel 
wide coverage of the knots. We therefore make an ex- 
trapolation of the results taken with these observations 
to the total area of the knots. 

The overwhelming majority of the PAH emission comes 
from the nucleus, which is also the region with the most 
emission at 8/im (Fig. From the star forming ring, 
Knots A, C, E, and F have PAH associated with the 
star clusters. Knot E is the one with the highest PAH 
emission, and this is also expected from the 8/iin IRAC 
image in Fig. \2§ ■ 

3.3. Ionic lines 

In Table [3] we present the fluxes of the main ionized gas 
lines in the low resolution spectra, as measured using 
SMART. The errors listed reflect the S/N ratio at the 
wavelengths where the lines were measured. With these 
fluxes we calculated the ionic line ratios [NeIII]/[NeII], 
[SIII]18.7/[SIII]33.5 fj,m, and also [Sill] /[SIII]. 

The [Nelll] / [Nell] rate, because it is a rate between 
two ionization states of the same element, it is sensitive 
to the effective temperature of the ionizing stars, and 
therefore can be used as a measure of the hardness of the 
radiation field. This ratio varies typically between 0.05 
and I in starburst galaxies and HII regions in normal 
gal axies, and is typically greater than 1 in dwarf galax- 
ies (|Thornlev et al.ll2000D . We could only measure the 
[NeIII]/[NeII] ratio for knots E, F, and for the nucleus, 
and the results are 0.52, 0.67, and 0.08. The ratios for 
knots E and F are compatible with very young massive 
clusters, and are similar to those found in very young 
clusters in other interacting systems such as the Anten- 
nae, where star forming regio ns with [Nelll] / [Nell] ratios 
between 0.30-0.73 are found ISniideri (|2007f) . The ratio 
for the nucleus is much smaller than the ratios for the 
ring knots and are similar to the ratios foun d in the nuclei 
of sta rburst galaxies, like NGC 253 (0.07; iDevost et al.1 
(l200l ). 

The [SIII] 18.7/ [SIII] 33 . 5 /mi ratio, as it is a ratio be- 
tween two lines of the same ionization state, but with 
different critical densities for collisions with electrons 
(|Osterbrockl H989). Therefore it is sensitive only to the 
density of the ionized ISM. The [SIII]33.5/jm line was 
detected in all knots, but the [SIII]18.7/im was detected 
only in knots E, F, G, and in the nucleus. As seen in 



Table H the values for [SIII] 18.7/[SIII]33.5 a™ ratio for 
these knots vary between 0.50 in knots E and F and 1.34 
in knot G although with a large measurement error for 
the last case. 

The [Sill] /[SIII] line r atio is used as a starburst /AGN 
diagnostic ([Dale et al.l l2006h . This is because the 
[Si II] 34.82 /im line is a significant coolant of X- 
ray ionized regions or dens e photodissociation regions 
(Hollenbach & Tielens 1999), more commonly associated 
with AGN, whereas the [S III]33.48 /jm line is a strong 
marker of H II regions. The typical value found in AGNs 
is around 3, whe reas for nuclear starbursts is around 1 
(|Dale et al.l 120061 ). We measured this ratio for all knots, 
varying from 0.82±0.22 in the nucleus to 1.97±0.51. The 
low signal-to-noise ratio in clusters B, C, and D is re- 
flected in the high errors of the [Sill] /[SIII] measure- 
ments, and therefore all the values for this ratio are com- 
patible with the values for nuclear starbursts found in 
([Dale et al.ll2006l ). 

3.4. High resolution spectrum of the nucleus 

To study the ISM conditions in the nucleus with 
greater accuracy we took a high resolution spectrum of 
the nucleus of NGC 2445, which is presented in Fig. [5] 
We scaled up the SH spectrum by multiplying it by a 
factor of 1.7 to match the LH spectrum flux at 20/im. 
This factor is much lower than the difference between 
the SH and LH slits, but that is simply due to the fact 
that the nucleus is practically a point source. The LH 
and SH spectra were then joined together, resulting in 
a single 10 — 35/im high resolution spectrum of the nu- 
cleus. The main features present are the ionic lines, like 
[Nell], [Nelll], [SIII], and [Sill]. PAH bands are also visi- 
ble, as the 11.3/imband, and the 16 — 1 8/im feature. Also 
present are the H2 rotational lines. These features allow 
a physical characterization of the ISM of the nucleus, and 
they can be used to diagnose electron temperature and 
density and H2 gas temperature. The measured line ra- 
tios are listed in Table [3J along with the lines measured 
in the same knot using the low resolution spectrum. We 
can see that the low resolution line fluxes have roughly 
half the flux as the high resolution lines. This is due to 
the difference in the extraction apertures. The area of 
the LH slit is 231 arcsec 2 , nearly twice the area of the 
low resolution extraction aperture chosen for the nucleus, 
which is 104 arcsec 2 . This gives a ratio of 2.22, whereas 
the ratio between low and high resolution apertures for 
both the [Nell] and [Nelll] flux is 2.04, meaning that 
the surface brightness of the lines decreases with a wider 
aperture. However, to do a comparison of the contin- 
uum fluxes one has to estimate the background flux. We 
measured the background flux from an off-source SL slit, 
which gives 116 mJy and 95 mJy for 15/im and 24/im, 
after scaling for aperture sizes. That means that the 
background accounts for ~ 55% of the high-resolution 
continuum at 15/im and ~ 32% of the high-resolution 
continuum at 24/im. The ratio of the continuum flux at 
15 /xm is thus 1.6, meaning that the surface brightness 
of the warm dust decreases with the aperture. The ratio 
for the 24/xm continuum is 2.2, meaning that the surface 
brightness of cooler dust does not vary with slit size. 

Comparing the ratios derived from high resolution with 
the low resolution ratios, we find a remarkable match for 
the [Nelll] / [Nell] ratio, being 0.08 for both spectra, but 
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for the [SIII] 18.7/ [SIII] 33.5 /im ratio there are significant 
differences, being 0.73 for the high resolution spectrum 
and 0.49 for the low resolution spectrum. This reflects 
the different flux ratios between high and low resolution 
apertures: [SIII]18.7/jm has a flux ratio of 2.38±0.16 thus 
maintaining its surface brightness, whereas [SIII]33.5/im 
has a flux ratio of 1.61±0.09 thus decreasing its surface 
brightness. This difference possibly reflects a decrease of 
the ionized gas density as the distance to the nucleus in- 
creases. Note that the LH slit could be contaminated by 
emission form the "shocked" gas region, and thus it can 
affect the line ratios that use lines [Sill] and [SIII]33.5/im 
and overestimate the nuclear H2 component. 

3.5. Photometry 

Photometry of each knot was performed using aper, 
an IDL tool from the IDL Astronomy Library. We used 
the FUV GALEX image to set an aperture size for each 
knot, as seen in Fig. [2 We used circular apertures of 
4". 5 for knots C, F, and the nucleus; 6" for knots B and 
D; and 7". 5 for knots A, F, and G. The sky background 
was determined using circular apertures of the same size 
as the ones used to measure each knot. With these aper- 
tures, we measured the flux within selected regions inside 
Arp 143 - but avoiding the knots - and averaged the re- 
sults. We applied extended source corrections for IRAC 
photometry. These corrections depend on the size of the 
aperture and do not exceed 6.8%. The resulting fluxes 
of each of the knots are listed in Table |U 

Arp 143 was not observed with the MIPS instrument 
at 24/im. However, in order to allow us to use some 
of the well known star- formation indicators (e. g. 
ICalzetti et al.l (|2007l ) ) associated with the knots based 
on MIPS observation, we decided to create pseudo-MIPS 
24/tm photometric points by extracting 24/im data from 
the full-IRS LL map (see Fig. [2]-last panel) using the 
CUBISM-extracted spectra, and the MIPS 24/im filter 
response curve. These values are tabulated in Tabled] 

4. RESULTS 

Our goal is to characterize the physical origin of the 
H2 line emission, its role on the appearance of the ring of 
star formation around the nucleus of NGC 2445, describe 
how the ISM in the knots evolves as they age. To achieve 
this goal we now analyze the observations described in 
the past chapter according to the following points: a) the 
study of the morphology of the system, describing it in 
the present state, in order to compare it to the present 
theories on the evolution of ring galaxies; b) the study of 
Hi excitation region and comparison with predictions, in 
order to study its physical origin; c) the ISM properties 
of the knots and its evolution, using the mid-infrared 
spectra to study the ISM in the knots, and UV to K-band 
photometry to study their ages; and d) the star formation 
rates in the knots, and the future of star formation in the 
system. 

4.1. Morphology 

In Fig. [7| we present a composite of FUV (blue), V- 
band (green), and 8.0 /xm emission (red), which traces 
emission from hot dust. Most of the luminosity in the 
V-band comes from the nuclei of the two galaxies, NGC 
2444, and NGC 2445, and a ring of star formation knots 
that surrounds the nucleus of NGC 2445. The nucleus 
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Fig. 7. — Composite of FUV (blue), V-band (green), and 8 iim 
(red) images of Arp 143. 

of NGC 2445 dominates the emission in the infrared, in- 
dicating the presence of large amounts of warm dust, 
but it also has a large unobscured population of main 
sequence and old stars, as it has also significant opti- 
cal and near infrared components. For a more quan- 
titative comparison, the 8/im flux from from the nu- 
cleus is about 5 times the combined 8/im flux from the 
knots. The knots are located alon g the outer edge of an 
HI crescent (|Appleton et al.lfl992f ). The knots are com- 
posed by several massive young clusters, visible in blue, 
and some are surrounded by PAH emission, visible in 
red. There are some mid-infrared "arms" connecting the 
nucleus to some of these regions. These arms are con- 
nected to the collisional nature of the ring - spokes are 
expected in gas tha t is collecting downstream of the ring 
(|Struck et al.l ll996). The only other galaxy with simi- 
lar feat ures reported in th e literature is the Cartwheel 
galaxy ([Struck et al.lll99(a ). 

NGC 2445 is about 20 kpc across. Some peculiarities 
reveal a connection between the massive star clusters and 
the dust, like the easternmost knot at about 12 kpc from 
the nucleus, identified in the FUV image in Fig. [5] as 
knot A. In Fig. [7] this FUV emitting knot seems to be 
surrounded by 8/im mid infrared emission (in red), as 
it is brighter around the knot and less bright inside a 
radius of ~ 1 kpc from the center of the knot. This 
can be interpreted as a shell of warm dust heated by the 
massive stars. Other knots have 8 /tm counterparts, like 
knot B, C, E, and F. For knots C and F, the 8/im emis- 
sion can be seen at the same site as the optical cluster, 
whereas in knot E, the 8/im is found somewhat offset 
to the SE. Since the 8/im IRAC band is dominated by 
PAH emission features, this means that the knots can 
be divided into two groups: those with associated PAH 
emission (A, C, E, and F), and "bare" knots, for which 
little or no PAH counterpart is observed (B, D, and G). 
This confirms what is observed in the spectra of Fig. [5] 
where PAH features are barely observed in knot D. The 
nucleus of NGC 2444 also shows substantial 8 /tm emis- 
sion, reflecting substantial massive star formation. The 
FUV emission comes mainly from knots A, E, and G, 
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showing emission from young massive stars, whereas the 
IR-bright knot F emits little UV. An even more strik- 
ing difference between mid-infrared and UV brightness, 
and Fig. [2] shows it clearly, is seen in the nucleus, the 
brightest region in 8/im, with very little UV emission. 

The companion lenticular galaxy NGC 2444 is promi- 
nent only in the optical and K-band, and is rather weak 
at longer wavelength. This supports the idea that much 
of the gas has been swept out of the galaxy in the past 
and it is dominated by an old inactive stellar population. 
It is possible that a faint "arm" extended from NGC 2444 
behind NGC 2445, but this is hard to prove without kine- 
matic data. 

4.2. The Shock Front 

Rotational emission lines of H2 may arise through three 
different mechanisms: UV excitation in PDRs surround- 
ing or adjacent to the HII regions; shocks that accelerate 
and modify the gas in a cloud, collisionally exciting the 
H2 molecules; and hard X-ray photons capable of pene- 
trating the molecular clouds and heating large ionizing 
columns of gas. Several rotational H2 lines were detected 
in the wavelength range of our observations and they are 
listed in Table [1] Given that these lines trace the molec- 
ular gas at different temperatures, studying the spatial 
distribution of these excitation lines may give us a clue 
about the main excitation mechanism. 

In Fig. [5] we present contour maps of the H 2 S(0) (left) 
and H2 S(l) emission (right) in Arp 143, overlayed on 
the IRAC 8/im emission color map. The H2 emission 
in Fig. [5] is concentrated in two clumps-from the nucleus 
and in a crescent-shaped ridge (including knot F) west of 
the nucleus corresponding to the Hi-crescent seen in VLA 
observations. The emission also is similar, but not identi- 
cal to CO (1-0) interferometry distribution Higdo n et all 
1997). These latter observations trace the distribution 
of cool H2 emission through collisional excitation of the 
CO molecule. The main difference between the Spitzer 
(warm) H 2 observations and the CO (cool H 2 tracer) 
is that the warm gas seems to better define a crescent- 
shape. Indeed it is precisely here that a shock wave 
would be expected to heat up the H2, supported by the 
fact that the S(l) map more clearly defines the crescent 
shape than the S(0) line, which is dominated by cooler 
excitation. This is compatible with the theoretical sce- 
nario i llustrated in Fig. [5J taken from [Applcto n et al.l 
(|1987af) . which represents the density wave formed after 
an off-center collision between two galaxies (the site of 
the collision is marked with an X), which expands out- 
wards from the center and also grows in length. It is also 
worth noting that the area chosen for the spectrum of the 
gas ring, demarcated in Fig. [^coincides with the area of 
the shocked H2 emission. Therefore we can consider this 
spectrum, in Fig. [4] as representative of the shocked H2 
region. 

The shape of the H2 emission in Fig.[5]is promising but 
by itself does not settle the case for either PDR or shock 
origin. A way to investigate the origin of the excitation 
of rotational transitions of the H2 molecule is to make 
a spatial comparison between the H2 emission and PAH 
emission. The H 2 emission-line strength should track 
the PAH emission-line and IR continuum strength in the 
case where the most common H2 excitation mechanism 
is UV excitation at the PDR interface with star forma- 



tion regions. However, we see a remarkable transition in 
Arp 143, from the nucleus, where PAH emission is quite 
strong, to the non-nuclear H2, which remains strong in 
the H2 lines, but is very weak in PAH (and continuum) 

e mission. 

iRoussel et al.l (|2007[ ) calculated the logarithm of the 
average ratio of the power emitted in the sum of the S(0) 
to S(2) transitions to the power emitted by the PAH fea- 
tures within the IRAC4 band, for HII nuclei and Seyferts 
from the SINGS survey. The average ratio for HII nuclei 
-2.19 ± 0.10 and for the Seyferts is -1.80 ± 0.34. We 
measured these ratios for the knots where H2 lines from 
S(0) to S(2) have been measured. We take 7.7/xm and 
8.6/im as the PAH features inside the IRAC4 band. For 
the nucleus of Arp 143, the ratio is -2.26, which puts 
it in the average of the HII nuclei. However, for knots 
E, F, and, G, the ratios are -1.22, -0.90, and -0.94, re- 
spectively. This puts them well above the averages for 
HII nuclei and also for Seyferts. For the gas ring re- 
gion, which largely coincides with the H2 crescent, this 
ratio rises up to -0.36. So the difference we observe in the 
H2/PAH ratio between the nucleus and the ring region in 
Arp 143 is a strong indication that the H 2 emission ridge 
seen in Fig. [5] is mostly due to shocks. X-ray emission is 
unlikely to be a strong contribution to the H2 heating, 
since the ratio of the [Sill] / [SIII] lines in Arp 143 (even 
in the nucleus) suggests an insignificant contribution to 
the excitation of [Sill] from X-rays. (X-ra ys can sig- 
nifican tly enhance this ratio-as discussed by Dal e et alJ 
(2006)). This strongly suggests shock-excitation within 
the crescent-shaped structure. A very strong H2 emis- 
sion in the absence of strong PAH emission is seen in the 
shock wave in Stephan's Quintet (Applcton et al. 2006) 
and in the spectrum of 3C326 (jOgle et al.ll2007[ L where 
strong shocks are implicated. The lack of a correspond- 
ingly large continuum is also one of the characteristics 
of the shocked H2 seen in the radio galaxy sample, and 
the Stephan's Quintet, as well as low-excitation emission 
lines. 

We can interpret the excitation diagrams in Fig. [6] also 
in the context of the different origins of the H2 emission. 
The temperatures derived in the excitation diagrams are 
compatible with either a PDR origin or a shock origin for 
the H2 excitation. The H2 temperature derived for the 
nucleus is higher than the temperature for the knots as- 
sociated to the H2 ridge, and to the H2 ridge as a whole, 
still does not reach the level required to rule out PDR 
origin, and the comparison with the PAH emission is cru- 
cial. However, these temperatures are just average gas 
temperatures, and it is possible the presence of a warmer 
gas component in the H2 shock region. The column den- 
sity of the gas is lower in the H2 ridge than in the knots, 
which is expected, since the warm H2 emission from the 
ridge is more diffuse, whereas the H2 emission peak co- 
incides with knot F. 

4.3. The Properties of Gas and PAHs in the Knots 

The knots in the ring of star formation are very young 
and are expected to have similar characteristics to young 
massive clusters discovered in starburst galaxies. We use 
our observations to study the components of the ISM in 
these knots, which consist of the molecular gas, ionized 
gas, and dust. 

The properties of the warm molecular g presented 
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Fig. 8. — Contour maps of the H2 S(0) (left) and H2 S(l) (right) emission in Arp 143, overlayed on an IRAC image. Contour levels 
are at 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5 , 2, and 3 MJy/sr for H 2 S(0) and 1, 1.4, 1.8, 2, 2.5, 3, 4, 6, and 8 MJy/sr for H 2 S(l). 




Fig. 9. — Schematic scenario of NGC 2445 where the knots are 
situated at the leading edge of the density wave. The contours are 
from th e off-center collision model E from Fig. 7 in Applet on et al.l 
( 1987a), representing a stage of the expansion process of a density 
wave. The arrows represent the expansion of the wave. Solid con- 
tours indicate elevation above the initial unperturbed disk values 
and dotted contours represent depressed levels. The nucleus is 
marked with a square. The site of the collision is marked with a 
cross. 

in Sect. 14. 2\ can be compared wi th the properties of the 
cold molecular gas as studied in (Hig don et al.lll997l ) us- 
ing CO observations. CO probes molecular gas up to a 
temperature of 100 K, but provides no information on 
the "warmer" gas which might be more directly linked 
to the source of activity. However, we can estimate what 
fraction of the total gas mass lies in higher temperatures. 
The fraction H2 (warm) to H2 (cool) is 4.8% for knots F 
and G, 2.6% for the nucleus, but 64% for knot E, assum- 
ing a CO-H2 conversion factor for the LMC. On average, 
the fraction of warm to coo l H2 amounts to 10%, w hich 
is nearly the same found bv iRigopoulou et "all (|2002f ) for 
a sample of starburst galaxies. 

Forbidden ionic lines can be used as indicators of 
the conditions of the ionized gas in the knots and the 
shock region. The ratio [NeIII]/[NeII] is used as a mea- 
sure of electron temperature of the ionized gas; the 



[SIII]18.7/[SIII]33.5 fim ratio is used as a measure of 
electron density. These ratios are compared to mod- 
eled ratios der i ved u sing photoionization diagnostics by 
ISnijders et all (|2007t ). assuming Salpeter IMF, stellar 
mass cutoffs M„ p = 100M o and M low = 0.02M Q and 
solar mctallicity. 

From the [SilI]18.7/[SIH]33.5 /im ratio we can esti- 
mate the density of the knots for which these lines are 
measured. This ratio is insensitive to the ionization pa- 
rameter Q and age. As seen in Fig. [10] (right), for any 
Q and age, knots E, F, G, and the nucleus have elec- 
tron densities in the order of 10 2 — 10 3 cm~ 3 . This is 
the average electron densi ty encountered in the center o f 
the starburst galaxy M82 (|Forster Schreiber et al.ll2001h . 
However Knot E, for example has a diameter of ~ 1 
kpc, about two times t he size of the central region of 
M82 (Bcir ao et al.ll2008h . and it is composed of individ- 
ual clusters that could have much higher densities. 

Fig. [TU] in the left shows the evolution of the modeled 
[Nelll] / [Nell] with age, for a given value of the ionization 
parameter. Due to difficulties in measuring the [Nell] line 
in the low resolution spectra, as it is merged with a PAH 
feature at 12.8£im, we could only measure the ratio for 
knots E, F, and for the nucleus. We can see that the mea- 
sured ratios correspond to knot ages of 1 - 4.5 Myr. The 
difference in the [Nelll] / [Nell] ratio between the nucleus 
(0.08) an the ring knots (0.5-0.7) means that population 
of recently formed massive stars in the nucleus is older 
than the young bursts in the ring, having an estimated 
age of 5 - 6 Myr. The fit to the low resolution spectrum 
of the nucleus gives an Ay of 0.02, which is too low to 
have any effect on the ratios, and therefore extinction 
does not affect the ratios significantly. 

PAH molecules with different physical characteris- 
tics produce different bands. The 6.2/im, 7.7/im, and 
the 8.6/mi bands are produced preferentially by ionized 
PAHs whereas the 11. 3 fan band is produced primarily by 
neutral PAH molecules. Therefore, the ratio 11.3/7.7/im 
may indicate the effect of PAH ionization. We calculate 
this ratio for all the knots. Given the uncertainties in the 
PAH flux measurements, we could only derive a reliable 
11.3/7.7/im ratio for knots F and H. We cannot see sig- 
nificant differences between the two ratios, at least not 
bigger than the uncertainties for knot F. 
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4.4. The Ages of the Knots 

With the fluxes listed in Table H we built SEDs of 
the main knots, in order to further constrain their phys- 
ical characteristics. To achieve this, we compare the 
SEDs with synthetic s pectra modeled with Starburst99 
iLeitherer fc Heckmanl ({1995) . There are several parame- 
ters we need to consider when fitting the model SEDs to 
our observations: 

• The mode of star formation. The star formation 
may occur continuously (continuous star formation 
model, or CSF model) or in a single, almost in- 
stantaneous, burst (instantaneous burst model, or 
ISB). We will adopt the ISB model for the star- 
forming knots, as there is no evidence that star 
formation has been ongoing for any considerable 
time in these knots. 

• The stellar initial mass function (IMF). We adopt 
a Kroupa IMF with an upper-mass cutoff of 100 
Mq and a lower mass cutoff of 0.02 Mq. 

• The initial gas metallicity (Z). We assume sub-solar 
metallicity Z = OAZp, , which is s i milar with the 
metallicity reported bv lJeske et al.l (|1986f h 

• The age of the star clusters. This is a free param- 
eter ranging from to 1 x 10 9 yr. 

• The effect of extinctio n. We use the Calze tti star- 
burst extinction law (|Calzetti et al.l [20001 ) on the 
synthetic spectra, 

• Dust emission. 

With these conditions, we fitted the synthetic spectra 
to the measured SEDs of each knot. The fluxes listed 
in Table |4] were converted to units of ergs _1 A -1 . We 
compensated the optical and UV fluxes for a galactic 
reddening o f E(B-V)=0.051, usi ng the extinction maps 
a nd laws of ISch lcgcl et al. (1998) and attenuation curve 
of lCardelli et al.l (|19 89) , the latter exclusively for the UV 
fluxes. Along with a pure stellar synthetic spectrum, the 
output of Starburst99 includes a nebular emission com- 
ponent as well, which was added to the stellar spectrum 
for knots C, E, and F. A set of free parameters, such 
as age, extinction, and continuum flux were adjusted for 
each knot in order to minimize the \ 2 value. The results 
are listed in Table El as well as the range of parameters 
within a 95 % confidence level. Knot B was left out, as 
knot B has a foreground star that interferes with the mea- 
surements. Although the nucleus could be modeled using 
the above assumptions, the measured optical to UV flux 
ratio appears always larger than the modeled spectrum. 
The variation between subsolar metallicity Z = O.4Z 
and solar metallicity does not modify significantly the 
parameter values. 

From the results of the SED fits in Table [5] we can di- 
vide the knots in two age groups: knots with ages from 
2-4 Myr (A, C, E, and F), and knots with ages from 
7 - 7.5 Myr (D, an d G). Given that the in teraction oc- 
curred 100 Myr ago lAppleton et all (|1987ah . this age dif- 
ference is not very significant in terms of the interaction 
dynamics. Probably due to a clumpy ISM, some clouds 
collapsed earlier than others, forming stars a little ahead 



in time than others. The estimated ages come into agree- 
ment with Fig. [10] (left), which predicted the same range 
of temperatures for knots E, and F. 

In Fig. [IT] we show two examples of SED fittings on 
knots C (with nebular emission) and D (without nebu- 
lar emission), being good representatives of the two age 
groups. We can see how the IRAC 5.8/im and 8.0/im 
fluxes, and the synthetic 24/im flux, deviate from the 
pure stellar model for knot D, indicating the presence of 
PAH and thermal dust emission. We can also see that 
the SEDs are not easily fitted with a simple instanta- 
neous starburst model, and an example is the excess of 
K-band emission in knot C, observed in Fig. [TTJ (left). A 
-ftf-band flux excess was also reported in previous stud- 
ies tha t fitted SEDs with Starburst99 (ISurace fc Sanderi 
(|1999h for warm ULIGs; iMazzarella et all (|1992l ) for the 
nuclei of Arp 220), which was compensated adding a hot 
dust (800 K) component contribu ting to 10-30% of the 
emission ([Surace fc Sandersl[l999l ). 

The nucleus of NGC 2445 was modeled using instanta- 
neous (ISF) and continuous star formation (CSF). Con- 
tinuous star formation is more a reasonable assumption, 
based on results from ring models, who predict a contin- 
uous fuelling of the nuclear starburst through gas inflow. 
In this case we used a Kroupa IMF with an upper-mass 
cutoff of 100 M Q and a lower mass cutoff of 1 Mq. We 
found that the best fit is achieved for the age of 500 
Myr. These ages are radically different than those ob- 
tained with [NeIII]/[NeII] ratios (5-6 Myr), and this 
is due to the fact that both [Nell] and [NcIII] lines are 
produced from the gas excitation by the youngest and 
hottest massive stars, and that ISF was assumed in the 
modelling of the lines. 

The age determination achieved here is a significant 
improvement in comp arison with past efforts, e. g. 
Applc ton et al.l (|1992f ). as their color-based determina- 
tions did not include UV bands. The FUV and NUV, as 
demonstrated here, are essential to the breaking of the 
age-extinction-metallicity degeneracy. 

4.5. Star Formation Rates 

With the H a , FUV and mid-infrared continuum fluxes 
at 15/mi, 24/im, and 30/im we can derive star forma- 
tion rates for each knot. Usin g the s tar formation rate 
calibration by iKennicutt et alJ (|1998h we can derive the 
star formation rate from the H a luminosity, which traces 
the very young and massive stars. However, the H Q 
flux is affected by internal extinction, and for the more 
dusty clusters, we sh ould use the Calzetti calibration law 
(jCalzetti et al.|[2~007f) , which includes the 24/im luminos- 
ity to compensate extinction that affects the H a emis- 
sion. Far ultraviolet is also a tracer of young massive 
stars, but even more affected by extinction than H Q . We 
use here the A v values in Table [5] to correct th e FUV 
fluxes and the calibration by ISalim et al.l (|2007f ) , which 
uses GALEX FUV calibrations. Mid-infr ared spectral 
diagno stics for the SFR were derived in iBrandl et al.l 
(2006), based on the fluxes at 15/im and 30/im. For 
each spectrum of the knots, we measured the flux at 14.5- 
15.5/m i and 29.5-30.5/im, a nd calculated the SFRs based 
on the IBrandl et al.l (|2006f ) calibrations. In Table [6] we 
list the corrected FUV, H a , 24/im luminosities, as well 
as the 15/im and 30/im fluxes. The SFRs derived using 
each one of these indicators are also listed. 
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Fig. 10. — Left: effect of the cluster age on the [NeIII]/[NeII] ratio. The model curves are computed for a cluster of 1O B M0, assuming 
Salpeter IMF, M up = IOOMq, and Mi ow = 0.2Mp). Each curve represents a different ionization parameter, and the solid curve is the one 
that approaches the value found in M82 by Forstcr Schreiber ct al. (2001), logU = —2.3. The horizontal lines indicate the [NeIII]/[NeII] 
for knots E and F, and the nucleus. The first dip in [NeIII]/[NeII] represents the ageing of the stellar population, after which WR 
stars are produced, increasing the [NcIII]/[NcII] ratio. The dip at 6 Myr occurs as most massive stars die through supernova explosions. 
Right: dependency of [SIII] 18.7/[SIII]33.5 /im ratio on the electron density. The curve is for a single ionization parameter approaching 
logC = —2.3. The horizontal lines indicate [SIII]18.7/[SIII]33.5 for knots E, F, and G, and the nucleus. Knots E, and F have the same 
value for [SIII]18.7/[SIII]33.5, which is 0.5. 



Adding all the SFRs of the knots, we arrive to a to- 
tal SFR of 2.16±0.17 Moyr^ 1 from H Q luminosities and 
2.40±0.35 M Q yr _1 using H Q + 24/zm lu minosities. This is 
similar to the 2.5 M Q yr _1 derived by (|Jeske et al.lll986l) 
using the H/3 fluxes but below the the far-infrared SFR of 
6.21 M^yr -1 calculated from IRAS colors ( jMoshir et alJ 
1990). It is also far below the SFR derived from the 
corrected FUV fluxes, which is in t otal 15.9 M^y r " 1 . 

The SFRs calculated by the iBrandl et all (|2006l ) 
method are systematically b elow the SFR calculated us- 
ing the (jCalzetti et al J 120071 ) calibrations, which in turn 
are below the SFRs derived from FUV. This can be ex- 
plained by the fact that most of the energy is detected 
in the optical and not absorbed by dust, which makes 
this system distinct from more archetypal starbursts like 
M82. Here, the UV and optical radiation from the clus- 
ters is heavily absorbed , making the A„ at least 10 times 
higher than in Arp 143 (|Beirao et al.ll2008[ ). This is miti- 
gated in case of the nucleus, which accounts by ~ 80% of 
the total 24fim flux, and thus the discrepancy is smaller. 
Notice that the star formation rate from UV is heavily 
dependent on the estimation of A„, which values within 
a 95 % confidence can vary by a factor of two, as seen in 
Tabled 

5. THE ROLE OF SHOCKS IN THE PROPAGATION OF 
STAR FORMATION 

In the past sections we discussed thoroughly the H2 
shocks occurring in NGC 2445 and the ages and den- 
sities of the knots that compose the star forming ring. 
The strong H2 emission reported this paper exists in and 
aro und the density wave, as traced by HI emission shown 
in (|Appleton et al.l 1992). This is expected from shocks 
created as the wave moves out through the disk. This 
is because the gas, unlike the stars in these models, can- 
not pass through each other and shocks should develop 
at caustics. In this section we will discuss the role of 
these shocks in the development of the knots, and its 
consequences on the morphology of Arp 143. 



5.1. The evolution of PAH emission with cluster age 

The SED modeling in this paper provided a further 
refinement of the age dating of t he knots, which is a m ain 
advantage over past attempts (|Appleton et al.lll992| ). It 
showed that not only the HII region nebulosity but also 
the clusters themselves are consistent with a recent very 
young burst. 

The evolution of PAH excitation with starburst age 
has be en speculated by other studies, like iRoussel et all 
(2005) in their analysis of dust excitation in NGC 300. 
We can try to observe this behavior taking advantage of 
our detailed diagnostics of knot ages. 

The improvement on the age dating of the knots al- 
lowed us to divide them into two groups: those with an 
age between 2-4 Myr old (knots A, C, E, and F), and 
those with an age of 7-8 Myr old (knots D and G). What 
is remarkable is that the younger group of knots cor- 
responds exactly to the knots where PAH counterparts 
are observed, whereas the older knots are those "bare" 
of PAHs, as discussed in Section 14.11 and seen in Fig. [7J 
This is the first time this ageing effect in star forming 
galaxies is observed in such chronological detail. 

The lack of PAH emission in older knots means that 
the PAHs have either "cleared out" , ceased to be excited, 
or destroyed during a timescale of 4-5 Myr. The first pos- 
sibility is that the PAH have been "swept" after ~ 6 Myr 
by winds from star formation regions. The 8/^m shell in 
knot A, for example, has ~ 3 kpc in diameter. Assum- 
ing a wind velocity of 300 km/s, the shell has been ex- 
panding for 5 Myr, a timescale comparable to the age of 
the knot. Another possibility could be PAH destruction. 
Small molecules could become more exposed over time 
to the hard UV field and are subsequently destroyed, as 
the dust shell expands and the PDRs become exposed to 
the outer radiation field. However, both these processes 
also occur in other galaxies where PAHs are observed at 
later ages, like the Antennae, so the cause for this be- 
havior resides probably on the nature of the ring galaxy 
system. Due to the shock wave traced by the H2 emis- 
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Fig. 11. — Spectral energy distributions for knots C (left) and D (right), with Starburst99 synthetic spectra. The error bars represent 
photometric measurement errors. 



sion, the knots in the ring were formed simultaneously, 
and therefore we assume an instantaneous burst for all 
the knots. This means that the UV radiation fades away 
after ~ 6 Myr, and no longer excites the dust surround- 
ing the knot. This can be observed in Fig. [TU1 where 
the [NeIII]/[NeII] ratio, a measure of the hardness of the 
radiation field, decreases dramatically after 6 Myr. The 
slope of the continuum of knots D and G longwards 25^m 
indicates the presence of cold dust, meaning that dust is 
present in these knots, but it is not heated by UV radia- 
tion. An older star-forming region with continuous burst 
like the nucleus of NGC 2445 keeps forming new massive 
stars which continue to excite PAHs and warm up dust. 
This could also be the origin of the SFR discrepancy seen 
in Sect. 14. 5| where the SFRs in the knots calculated from 
the infrared continuum are systematically lower than the 
SFRs calculated from H Q +24/im luminosities, but not 
the SFR in the nucleus. 

5.2. The H2 emission front and the simultaneity of knot 

formation 

The strong H2 emission reported this paper traces a 
density wave expanding outwards at a constant speed. 
In response to the density wave, the knots were formed 
simultaneously in situ, and that is confirmed by the nar- 
row age range found in Table 6. The fact that the H2 
follows the shape of the HI overdensity is further confir- 
mation that this is a coherent structure. 

If w e assume the shock velocity from (jHigdon et al.l 
Il997fl . 118±30 km s _1 and the distance of the clusters to 
the nucleus, we can calculate the propagation timescale 
of the shock from the nucleus outwards. Given that the 
knots are on an average distance of 10 kpc, the shock 
wave has been propagating from the nucleus since 85 
Myr ago. However, star formation has been occurring 
in the nucleus way before then, at least since 300 Myr 
ago. This is consistent with the ki nematic picture first 
advanced in lAppleton et al.l (|1992f ): a first interaction 
with NGC 2444, provoking the onset of star formation in 
the nucleus 300 Myr ag o, and the HI plume observed in 
lAppleton et al.l (| 198 7a): and a second interaction ~ 85 
Myr ago, provoking the emergence of a s hock wave which 
create s a gas overdensity visible in HI App leton et al.1 
(1992): lHigdon et al.l (|1997f) . leading ultimately to a ring 



of young star forming knots. These knots were created 
almost simultaneously, and as we have seen, blows out 
surrounding PAHs in a very short timescale, 5 — 6 Myr. 

Models of ring galaxies have shown the emergence of 
simultaneous star formation due to an exp anding shock 
wave a s a consequenc e of a head-on coll i sionlGerber et al.l 
(19961); iStruckl (11991 : lLamb fc HeaTnl ((2000). However, 
the timescales and details such as gas distribution, fre- 
quency of starbursts, etc., are strongly dependent on the 
position and kinematics of the galaxies, and there are no 
models for the particular case of Arp 143. 

The future of star formation in the ring of knots will de- 
pend on the gas reservoir. Given the mass of atomic and 
molecular gas, and the SFR for all the knots, we can cal- 
culate the dur ation of star fo r matio n in Arp 143 assuming 
constant SFR. lHigdon et al.l (|1997| ) found 1.25± x 10 9 M Q 
of atomic gas in Arp 143, and 2.2 x 10 9 M Q of cold molecu- 
lar gas. Taking the total gas content and an average SFR 
of 2.16 M Q , we calculate that the star formation in Arp 
143 will last for ~ 2 Gyr. However, the evolution of Arp 
143 as a ring galaxy is also likely to increase its SFR, 
meaning that this is an upper limit for the duration of 
star formation. 

These estimates do not, however, take into account the 
dynamical evolution of the system. The intruder galaxy, 
NGC 2444, is likely to swing back, and disrupt once more 
NGC 2445. A dynamical model will be extremely useful 
to predict the future of star formation in this system, 
and therefore in similar systems at higher redshifts. 

6. CONCLUSIONS 

Mid-infrared observations of Arp 143 have been pre- 
sented in this paper, along with ancillary data including 
GALEX UV images. Multiwavelength images from UV 
to mid-infrared show a bright dusty nucleus surrounded 
by young star forming knots. The four main conclusions 
of this study are the following: 

• Spectral line maps of the H2 rotational lines show 
a ridge of warm H2 emission that curves between 
the nucleus and the western knots. The H2 line 
flux related to PAH flux is nearly 10 times higher 
in the ridge than in the nucleus. The flux ratios 
between the sum of H2 S(0) and S(2) lines over the 
sum of PAH 7.7/im and 8.6/im lines reveal that this 
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H 2 ridge observed in Fig. [8] arises from shocked gas 
behind the wave that provoked the onset of the ring 
of star forming knots. This is one of the few cases 
were this kind of feature is seen. 

• With the use of photometry and fitting the SEDs, 
we improved greatly the age determination of the 
knots. The knots in the ring are all very young, 
varying from 2-7.5 Myr old and the nucleus is about 
500 Myr old. The ring of knots are a product of a 
shock wave that has been expanding from the nu- 
cleus since ~ 85 Myr ago, and were formed simul- 
taneously in situ. Behind the shock the molecular 
gas condenses, and it is traced by the H2 emission 
ridge. The improvement of age determination is 
achieved in this study especially due to the GALEX 
FUV and NUV bands, which are crucial to break 
the age-extinction-metallicity degeneracy. 

• The distribution of ages of the knots correlate with 
the presence of PAH emission. Younger 2-4 Myr 
old knots are associated with PAH emission shells, 
whereas older 7-8 Myr knots contain little or no 
PAH emission. The most plausible explanation for 



this would be an effect of the ageing of the massive 
stars created in a single instantaneous burst all over 
the ring, within a timescale of ~ 6 Myr, after which 
the UV radiation is no longer energetic enough to 
excite the PAHs. 

• Given the current reservoir of molecular gas, and 
assuming that the current star formation rate 
maintains itself constant, the star formation in Arp 
143 will last for about 2 Gyr. However, as it is very 
likely that the SFR will increase dramatically, this 
can only be an upper limit for the duration of the 
starburst. 



We would like to thank Y. Mayya for providing the 
H Q image. We also thank B. Groves and L. Snijders 
for making the Starburst99 and Mappings models avail- 
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Spitzer Space Telescope, which is operated by the Jet 
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TABLE 1 

Molecular Hydrogen Line Fluxes (Top Row), Equivalent Widths (Bottom Row) and Diagnostics of 

the Main Knots 3, 



Knot 


H 2 S(3) 


H 2 S(2) 


H 2 S(l) 


H 2 S(0) 


T (K) b 


Nff 2 c 


Mff 2 d 


L(H 2 )° 


A 






0.27 


0.27 










B 


0.22 


0.31 


0.39 


0.18 










C 


0.33 


0.38 


0.34 


0.50±0.14 


















0.604 










D 


0.39 


0.38 


0.27 


0.20 










E 


0.41 


0.31 


2.48±0.03 


1.12±0.08 


137 


30.6 


10.8 


1.35 








0.335 


0.340 










F 


1.85±0.10 


1.17±0.21 


4.24±0.32 


0.88±0.12 


280 


5.8 


0.92 


3.11 




0.386 


0.183 


1.106 


0.194 










G 


0.34 


0.32 


4.27±0.52 


0.70±0.12 


178 


10.4 


3.70 


1.78 




0.355 




1.285 


0.233 










Ring 


3.37±0.34 


1.57±0.18 


14.5±0.6 


5.70±0.69 


236 


3.1 


4.45 


7.46 




0.114 


0.027 


0.309 


0.172 










Nucleus 


1.66±0.05 


1.20 


3.60±0.15 


0.91 


309 


2.67 


0.42 


2.13 




0.022 




0.030 












Nucleus (high resolution) 




2.72±0.28 


5.94±0.36 


0.87 


330 


1.58 


0.60 


3.98 






0.008 


0.019 













a Fluxes in units of 10 21 W cm , and the equivalent widths, below the fluxes, are in fim. Fluxes without errors are 
upper limits. b Temperatures based on the following slopes: S(0)-S(1) for knot E, S(0)-S(3) for knot F, S(0)-S(1) for knot 
G, and S(l)-S(3) for the nucleus Column densities in units of 10 19 cm" 2d Warm molecular gas mass in units of 10 7 
Mq c Luminosities in units of 10 40 erg s _1 



TABLE 2 

PAH FLuxES a (top row) and Equivalent Widths (bottom row) in the Main Knots 



Knot 


6.2 (im 


7.7 fim 


8.6 fim 


11.3 (im 


17 fim complex 


11.3/7.7 


B 


1.186±0.395 


1.717±0.466 


0.574±0.396 


0.567±0.256 


0.146±0.103 


0.33±0.33 




1.14 


2.49 


0.949 


1.20 


0.045 




C 


0.878±0.229 


1.952±0.641 


0.371±0.236 


0.627±0.128 




0.32±0.26 




1.38 


3.02 


0.652 


1.73 






D 


0.640±0.359 


1.397±0.819 




0.446±0.392 




0.32±0.32 




6.49 


11.1 




2.93 






E 


2.083±0.814 


5.824±2.675 


1.256±0.805 


1.503±0.455 


0.665±0.632 


0.26±0.26 




4.42 


9.03 


1.72 


1.74 


0.756 




F 


0.922±0.304 


4.369±0.724 


0.734±0.334 


0.972±0.184 


0.635±0.228 


0.22±0.10 




2.93 


9.08 


1.26 


1.51 


1.17 




G 


0.683±0.580 


3.98±1.449 


0.724±0.416 


0.697±0.631 




0.18±0.18 




7.89 


20.7 


1.46 


1.25 






Nucleus 


22.65±0.441 


89.52±0.865 


13.39±0.406 


15.93±0.312 


7.346±0.353 


0.18±0.00 




3.04 


14.8 


2.12 


1.87 


0.595 




Ring 


1.14±0.36 


4.23±1.62 


0.74±0.34 


1.04±0.18 


0.89±0.31 






1.87 


7.97 


1.38 


1.60 


0.94 





a Fluxes in units of 10~ 20 W cm~ 2 

Struck, C, Appleton, P. N., Borne, K. D., & Lucas, R. A. 1996, 
AJ, 112, 1868 

Struck-Marcell, C, Scalo, M., & Appleton, P. N. 1987, Star 

Formation in Galaxies, cd. C. J. Lonsdale (NASA Pub. 

2466)(Washington:GPO), 435 
Surace, J. A., & Sanders, D. B. 1999, ApJ, 512, 162 
Theys, J. C, & Spiegel, E. A. 1977, ApJ, 212, 616 
Thornley, M. D., Forster Schreiber, N. M., Lutz, D., Genzel, R., 

Spoon, H. W. W., Kunze, D., & Sternberg, A. 2000, ApJ, 539, 

641 



Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 

Toomre, A. 1978, IAU Symp. 79, The Large Scale Structure of the 

Universe, ed. M. S. Longair & J. Einasto(Dordrecht:Reidel), 109 
Verontsov-Vclaminov, B. A. 1959, Atlas and Catalogue of 

Interacting Galaxies Part I(Moscow:Moscow Univ. Press) 
Verontsov-Velaminov, B. A. 1977, A&A, 28, 1 
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White, S. D. M. 1979, MNRAS, 189, 831 



TABLE 3 

Ionic Line Fluxes and Ratios of the Main Knots 



Knot 


fSIVl 

[Ol V J 


[Nell] 


[Nelll] 


[SIII]18/im 


[SIII]33/im 


fSilll 


[NcIII]/[NcII] 


[SIII]18.7//m/[SIII]33.5/im 


[SiII]/[SIII] 


A 






0.53 




0.99±0.09 


0.87±0.05 






0.87±0.15 


B 


0.24 




0.29 


0.43 


0.59±0.10 


1.16±0.06 






1.97±0.51 


C 


0.36 




0.33 


0.36 


0.43±0.12 


0.87±0.11 






1.14±0.90 


E 


0.90±0.04 


2.25±0.09 


1.16±0.04 


1.24±0.27 


2.48±0.29 


2.88±0.16 


0.52±0.04 


0.50±0.19 


1.16±0.23 


F 


0.34±0.02 


0.92±0.05 


0.62±0.06 


0.85±0.19 


1.70±0.22 


1.64±0.08 


0.67±0.11 


0.50±0.20 


0.96±0.20 


G 


0.84 


0.84 


0.90 


1.85±0.35 


1.38±0.21 


2.05±0.25 




1.34±0.63 


1.49±0.48 


Nucleus 


0.06 


30.8±1.36 


2.54±0.08 


8.71±0.25 


14.2±1.2 


14.56±1.60 


0.08±0.01 


0.61±0.08 


0.82±0.22 


Nucleus (high resolution) 




63.6±1.6 


5.17±0.60 


20.7±0.8 


28.4±5.3 


31.9±1.6 


0.08±0.01 


0.73±0.20 


1.12±0.21 


Ring 


1.12±0.14 


3.29±0.38 


0.99 


5.21±0.89 


10.0±1.12 


8.10±0.68 




0.64±0.15 


0.81±0.14 



TABLE 4 

Broadband Fluxes of the Main Knots 



02 



Knot R.A. (7 46) Dccl. (+39) 



FUV a 



NUV a 



ya 



R a 



K a 



3.6 /^m a 



4.5 (im a 



5.8 (im a 



1.0 /im a 



24 (im a 



A 
B 
C 
D 
E 
F 
G 

Nucleus 



58.21s 
55.38s 
54.22s 
53.32s 
52.94s 
53.97s 
54.22s 
55.12s 



00 min 51.24 sec 
00 min 25.75 sec 
00 min 21.25 sec 
00 min 18.25 sec 

00 min 51.25 sec 

01 min 06.25 sec 
01 min 18.25 sec 
00 min 55.75 sec 



0.257±0.005 
0.171±0.003 
0.093±0.002 
0.197±0.003 
0.348±0.005 
0.069±0.002 
0.167±0.004 
0.049±0.005 



0.311±0.019 
0.231±0.012 
0.113±0.008 
0.208±0.013 
0.409±0.019 
0.077±0.007 
0.224±0.018 
0.132±0.025 



0.504±0.032 

3.17±0.07 
0.323±0.026 
0.632±0.040 
0.754±0.050 
0.120±0.010 
0.636±0.047 
1.41±0.053 



0.544±0.035 

3.62±0.08 
0.360±0.029 
0.755±0.041 
0.895±0.042 
0.143±0.012 
0.808±0.040 
2.48±0.064 



0.587±0.031 
4.00±0.08 
0.444±0.028 
0.916±0.043 

1.08±0.04 
0.173±0.020 
0.908±0.037 
3.85±0.073 



0.761±0.03 
6.18±0.06 

0.644±0.02 
1.25±0.03 
1.39±0.04 

0.309±0.02 
1.25±0.03 
16.5±0.1 



0.400±0.172 

2.74±0.32 
0.501±0.274 
0.587±0.200 
0.998±0.246 
0.318±0.142 
0.560± 
10.7±0.6 



0.257±0.186 

1.73±0.29 
0.308±0.139 
0.276±0.173 
0.647±0.337 
0.193±0.122 
0.335±0.129 
7.69±0.52 



0.878±0.447 
2.03±0.43 
1.14±0.32 

0.604±0.368 
2.78±0.51 
1.10±0.32 

0.484±0.279 
29.5±1.0 



1.97±0.74 
2.69±0.65 
2.71±0.54 

0.645±0.597 
6.51±0.84 
2.91±0.55 

0.841±0.719 
87.6±1.8 



Fluxes in units of mjy Errors represent Ict deviations from the median flux between 22 — 28/^m. 
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TABLE 5 

Properties of Knots Derived from Starburst99 
with Best Fit (Upper Row) and 95 % Confidence 
(Bottom Row) 



Knot 


Age (Myr) 


A v 


Mass (10 6 M©) 


A 


3.5 


1.0 


12.6 




O A 

Z-4: 


0.8-1.2 


10-15.8 




O.O 


1.0 


1 O ft 




1-5.5 


1.1-1.9 


7.9-15.8 


D 


7.5 


0.5 


12.6 




6.5-8.5 


0.4-0.6 


10-15.8 


E 


3.5 


1.2 


25.1 




3.5-5 


1.1-1.4 


20-25.1 


F 


2.5 


1.3 


6.3 




0.5-6.5 


0.8-1.8 


4.0-10 


G 


7.5 


0.6 


15.8 




6.5-8.5 


0.3-0.8 


12.6-20.0 


Nucleus (ISF) 


16 


3.0 


1122 


15-16 


2.6-3.4 


1120-1260 


Nucleus (CSF) 


500 


2.7 


6.3 




400-800 


2.2-3.1 


5-7.9 



TABLE 6 

Star Formation Rates of the Main Knots 



Knot L(FUV) a L(H Q ) b L(24^tm) b F(15Atm) c F(30Atm) c SFR(FUV) d SFR(H a ) d SFR(H a +24^m) d SFR(F(15)+F(30)) d 



A 1.63±0.10 14.3±0.2 237±11 2.8±1.3 6.7±1.0 1.76±0.06 0.11±0.02 O.lliO.Ol 0.05 

B 8.97±1.83 127±6 1.9±0.7 4.0±0.7 0.07±0.02 0.07±0.01 0.03 

C 1.87±0.04 16.7±2.5 107±5 1.5±0.4 2.7±0.7 2.02±0.04 0.13±0.02 O.lliO.01 0.02 

D 0.390±0.006 1.1±0.1 1.6±0.1 0.419±0.06 0.012 

E 3.49±0.05 45.2±4.1 491±10 4.4±1.1 15.0±1.6 3.77±0.05 0.36±0.03 0.32±0.02 0.10 

F 0.880±0.025 27.6±3.2 365±9 3.5±0.9 15.0±0.8 0.951±0.027 0.22±0.02 0.21±0.01 0.10 

G 0.421±0.010 12.8±2.2 267±10 2.6±1.9 5.6±1.2 0.463±0.011 0.10±0.02 O.lliO.01 0.04 

Nucleus 6.04±0.62 67.9±4.96 8420±363 62.9±1.8 274.9±4.8 6.53±0.67 0.54±0.04 1.46±0.27 1.63 



a Luminosities in units of 10 28 erg s 1 Hz b Luminosities in units of 10 39 erg s lc Fluxes in mjy d Star formation rates in Mq yr 



